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ABSTRACT 

We study the physics behind the neutrino-driven mechanism for gamma-ray bursts and hy- 
pernovae, deriving the critical density at which these outbursts occur in the collapsar model. 
The agreement between this derivation and results from past collapsar simulations (MacFadyen 
& Woosley 2000) is excellent, implying that we have captured the essential physics. We then 
use this derivation to study a range of progenitors for collapsar gamma-ray bursts. We derive 
how much of the star will accrete onto the black hole core before the infall density drops below 
this critical density, leading to an estimate of the remnant black hole mass for GRBs and hyper- 
novae. We also estimate the time delays between gravity wave or neutrino signals and the onset 
of the explosion or burst event. This derivation, combined with future observational constraints, 
provides a physical insight into the structure of the GRB progenitor. 



Subject headings: stars:evolution 



Introduction 



The leading model for long gamma-ray bursts 
(GRBs) and hypernovae 1 assumes that they are 
driven by energy released from a disk accreting 
onto a black hole resulting from a massive stel- 
lar collapse. This accretion energy is converted to 
the explosion energy in the GRB/hypernova jet 
either through some magnetic field mechanism or 
through the annihilation of neutrinos emitted by 
the hot disk. In this paper, we study the physics 
behind the neutrino driven explosion mechanism, 
deriving critical densities at which an explosion 
is launched. For long-duration (> 10s) GRBs 
and all hypernovae, the progenitor is a massive 
star that either doesn't or weakly explodes via 
the normal supernova mechanism: a.k.a. a col- 
lapsar (MacFadyen & Woosley 2000, MacFadyen, 
Heger, & Woosley 2001). By studying the struc- 



1 We restrict the definition of hypernova in this paper to 
those jet-driven explosions (see Heger et al. 2002) that are 
not associated with gamma-ray bursts. 



ture of these massive stars at collapse and com- 
paring to our critical densities, we can estimate 
the delay between stellar collapse and the launch 
of the jet as well as the ultimate black hole rem- 
nant mass. Future observations of these delays 
and remnant masses provide constraints on the 
progenitor and/or jet mechanism for gamma-ray 
bursts. 

2. Neutrino-driven Explosions 

Neutrinos emitted from black hole accretion 
disks push out against the infalling star both 
through the direct absorption and scattering of 
neutrinos on infalling matter and through the 
energy deposited in it as electron/positron pairs 
and photons from neutrino annihilation. Neu- 
trino annihilation, which provides most of the ex- 
plosion energy, occurs primarily above the black 
hole along the rotation axis of the accretion disk. 
Popham et al. (1999) studied a range of likely ac- 
cretion disk structures and calculated their resul- 
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tant energy deposition through neutrino annihila- 
tion. Using this study, we can derive the value of 
the matter density above the black hole for which 
the energy deposited by neutrinos is able to drive 
an explosion, and obtain an estimate of the mass 
of the black hole remnant left behind by this ex- 
plosion. 

We will define the onset of the explosion as the 
moment when neutrino momentum deposition ex- 
ceeds the pull of gravity on matter flowing onto 
the accretion disk along the rotation axis. The 
scattering opacity (k sc ) for neutrinos is roughly 
(Janka 2001): 

^« 5a ! + 1<r ? < l > -g-(y, + y P ), (i) 
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where m u w 1.66 x 10~ 24 g is atomic mass unit, 
m c c 2 = 0.511MeV is electron rest-mass energy, 
(To = 1-76 x 10~ 44 cm 2 , e„ is the neutrino energy, p 
is the density above the rotation axis, Y n — n n /nb 
and F p = n p /n^ are the number fractions of free 
neutrons and protons. The corresponding absorp- 
tion opacity (/c a b) for neutrinos is (Janka 2001) 



3a 2 + 1 cr < (-1 > P 



(m e c 2 ) 2 



(Yn,Y p ). (2) 



For the conditions in our disks, electron capture 
produces over half of the total neutrinos and elec- 
tron neutrinos arc the most abundant neutrino 
species. Assuming this species dominates the ab- 
sorption and scattering, the total neutrino opacity 
is 

Kt w 1.5 x 10- 17 p(^ e / 4MeV ) 2cm_1 - ( 3 ) 

Here we have used the following assumptions: 
axial-vector couplings set to the charged-current 
axial-vector coupling constant in a vacuum, 
a = —1.26, Y n w Yp w 0.5, the neutrino tem- 
perature (T„ c ) is related to neutrino energy by 
21(fcT„ e ) 2 =< e 2 > and that the entire neutrino 
flux is in the electron neutrino species. Had we as- 
sumed an equal mix of electron and anti-electron 
neutrinos, the total opacity would not be differ- 
ent by more than 20%. As a shell of matter falls 
towards the black hole, it is supported by the mo- 
mentum of the scattered and absorbed neutrinos: 



K t dr L v 1.5 x 10- 17 (kT„j4MeV) 2 L„ 
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where L v is the neutrino luminosity, m s hcii = 
4pnr 2 dr is the mass, dr the thickness, r the radius 
of the shell, c is neutrino velocity w the speed of 
light. 

The corresponding acceleration from neutrino 
annihilation requires detailed models of the black 
hole accretion disk system. The disk models 
and annihilation calculations from Popham et 
al. (1999) led to a value for the energy de- 
posited along a surface per unit path length 
([e] = crgss -1 cm -1 ). With this deposition rate 
e(r), which is a function of height above the disk, 
we can calculate the acceleration due to neutrino 
annihilation inside a 30° cone along the rotation 
axis: 

e{r)dr 1 e(r) 
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The onset of the explosion occurs when a g + 
a v + QW > where a g — — GM^wjr 2 is the grav- 
itational acceleration with gravitational constant 
G, and black hole mass Mbh 2 - Using equations 4 
and 5, the acceleration condition for the onset of 
the explosion can be translated into a threshold 
condition on the density: 



Pcrit < 



4e(r) 



47tGMb H c- 1.5 x 10- 17 (fcT„ /4MeV) 2 L„' 

(6) 

For our estimates, we take the radii where the en- 
ergy deposition from neutrino annihilation peaks 
(~20km). Figure 1 gives critical densities for 
set values of e versus an effective luminosity 
[(fcT I , o /4MeV) 2 L l/ ]. The critical density rises 
sharply as the "Eddington" luminosity for neu- 
trinos is reached. The data (e at 20 km and L v ) 
for a scries of disks from Popham ct al. (1999) are 
also shown for comparison. 

These critical densities correspond to an accre- 
tion rate inside the 30° cones along the rotation 
axes of 

M cr:i t = 0.5367rr 2 p crit Wff (7) 

where r — 20 km and vg = ^/2GMbhJt. The ac- 
cretion rates for massive stars are initially much 
higher than these critical values, but the accretion 
rate decreases as material from increasingly out- 
ward layers of the star falls onto the black hole. In 
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2 Note that we have neglected the disk mass which, in 
Popham ct al. (1999) can be as high as 0.5 Mq. Even 
for this high disk mass, the change in the gravitational ac- 
celeration is less than 20%. 
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figure 2 we plot the critical accretion rates for disks 
accreting 0.1 and 10 Mqs^ 1 . Note that the accre- 
tion rate through the disk need not be related to 
(and at late times, is much higher than) the accre- 
tion along the rotation axis. For comparison, we 
plot the accretion rates of stars with masses rang- 
ing from 20-60 M (Rauscher et al. 2002) versus 
mass coordinate of the accreting mass layer. 

When the accretion rate of stellar material 
falling along the polar axis drops below the critical 
accretion rate for a given disk, the neutrinos from 
the disk will drive an explosion and, very likely, 
disrupt the entire star. For a given disk structure, 
then, this crossing point gives a rough estimate of 
the remnant mass (e.g. from Fig. 2 we see that 
for a disk accretion rate of 0.1 M Q s~ 1 onto a black 
hole rotating near breakup: a=0.95, the remnant 
black hole mass for a 3OM star is roughly 10M Q .) 
Note also, that if the neutrino emission was as high 
as that given by the disks accreting at 10M Q s _1 
from Popham et al. (1999), then the explosion 
would occur immediately. However, Di Matteo, 
Perna, & Narayan (2002) have found that these 
disks actually produce a much lower neutrino flux. 
Even if such high neutrino luminosities could be 
constructed, they would disrupt their star imme- 
diately. 

We note that these results are only rough ap- 
proximations, since a number of caveats limit the 
quantitative accuracy of these calculations. Be- 
sides assuming a fairly simple model for the mo- 
mentum deposition, relativistic effects on the neu- 
trino energy and momentum, the effects of rota- 
tion on the infalling matter and the evolution of 
the black hole mass and spin were all neglected. If 
magnetic fields drive the explosion, the explosion 
can occur at much earlier times, and similarly, the 
disk winds seen by MacFadyen & Woosley (2000) 
also would eject the star at earlier times. Most of 
these effects would lead to higher critical densities 
and lower remnant black hole masses. However, if 
we restrict ourselves to neutrino-driven explosions 
only, the major uncertainty is the structure of the 
progenitor star. Figure 3 shows the infall accre- 
tion rates for rotating and non-rotating stars of 
40 and 60 M (Heger 2002). Note that the rem- 
nant mass differs dramatically depending upon the 
amount of rotation. Comparing the non-rotating 
60 M Q case of Fig. 3 with that of Fig. 2 shows the 
variations arising from stars produced with differ- 



ent versions of the same stellar code: Rauscher 
et al. (2002) and Heger (2002) use different ver- 
sions of the Kepler code (Weaver, Zimmerman, & 
Woosley 1978). In addition to uncertainties in sin- 
gle star evolution, it is likely that collapsar GRBs 
arise from binary systems, and binary progenitors 
of GRBs have not yet been constructed. 

3. Conclusions 

For most values of the disk accretion rates 
(within roughly 0.05 — lM Q s _1 ) the critical den- 
sity lies in a fairly narrow range between 10 4 — 
10 8 gcm -3 . These densities agree well with the 
results from MacFadyen & Woosley (1999). For 
lower disk accretion rates, the energy deposition 
from neutrino annihilation decreases dramatically 
and the critical density for explosions drops be- 
low 10gcm~ 3 . The material along the poles does 
not reach this density until nearly all of the star 
has accreted onto the black hole, and this will 
not produce a GRB. Similarly, higher neutrino 
luminosities would disrupt the star immediately 
and probably never form in nature. With such a 
narrow range of explosion conditions at the black 
hole source, one would expect some similarities in 
the GRB outbursts produced by neutrinos. How- 
ever, bear in mind that the propogation of the 
jet through the star (which can differ in different 
progenitors) can significantly alter the observable 
outburst. 

What does differ dramatically is the time after 
the collapse at which the explosion occurs. From 
figures 2 and 3, we know the mass zone of the 
star that is driven to explosion by the neutrinos. 
Using the 40, 60M© of Heger (2002) 3 and assum- 
ing that, along the rotation axis, the critical mass 
zone collapses at free-fall, we find that disks ac- 
creting at 0.1M Q s _1 drive explosions at much dif- 
ferent times. For the 40M Q star, infalling material 
along the pole is turned around at 50,60 or 10 4 s 
after the initial collapse for black hole spin rates 
of a=0.95,0.75 and 0, respectively. For the 60M Q 
star, the corresponding explosion times are: 35, 
300, 3.1 x 10 6 s. Hence, we should expect a consid- 
erable delay (at least ~ 30 s) between the collapse 
(as signalled by the initial neutrino and gravita- 
tional wave signal) and the launch of the GRB 

3 These rotating stars are the closest progenitors for GRBs 
and hypernovae yet constructed. 
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explosion. 

There are, however, upper limits to the delay. 
The explosion must occur, after all, before the 
disk accretes. Using the convection dominated ac- 
cretion flow (CDAF) solutions from Narayan, Pi- 
ran, & Kumar (2001), we find that for accretion 
rates above 0.05M Q s -1 and a total mass accreted 
through the disk less than 10M Q , the total time in- 
cluding the delay since the collapse plus the burst 
duration can not be higher than 200s. Those disks 
which do not explode before this time will not ex- 
plode before the disk dissipates. 

Requiring that the accretion rate is at least 
O.O5M0S -1 places an upper limit on the disk for- 
mation radius (the radius at which the angular 
momentum in the infalling stellar material sup- 
ports it against the gravitational pull of the black 
hole) . Similarly, requiring that the disk maintains 
its accretion rate long enough to survive the de- 
lay times given above places a lower limit on the 
disk formation radius. Using the CDAF solutions 
from Narayan et al. (2001) and assuming that the 
disk a viscosity ~ 0.01 -0.1, M D > 0.05MQS" 1 , 
^accretion > 40s, a total mass accreted through the 
disk = 10M Q with at least 1M Q in the disk, yields 
a range of disk formation radii from 28-7000 km. 
This corresponds to specific stellar angular mo- 
menta in the range: 3 x 10 16 — 5 x 10 17 cm 2 s _1 . 
Although this is a narrow range of stellar angular 
momenta, it does lie within many of the current 
rotating models (e.g. Heger 2002). 

For neutrino-driven explosions from black hole 
accretion disks, we can derive the remnant mass 
of black holes. Stars between ~ 20 — 40M© are 
likely to have weak supernova explosions which 
ultimately lead to considerable fallback and the 
formation of a black hole (Fryer 1999). If these 
stars have insufficient angular momentum, they 
produce a range of black hole masses between 
2 - 15 M0 (Fryer & Kalogera 2001). If instead, 
they are rapidly rotating, they can form hyper- 
novae and lower mass black holes (Nakamura et 
al. 2000). Simulations by MacFadyen, Woosley, 
& Heger (2001) show that ~ 10 4 s after their weak 
explosions, the infall rates of these stars will drop 
below our critical densities (assuming disk accre- 
tion rates in the range: 0.05 — \Mqs~ 1 ) and will 
drive explosions. These "collapsar type II" objects 
will range from 2 — 5M©. 

For direct collapse black holes which are the 



more likely GRB candidate, the currently most 
reliable progenitors (rotating 40,60 M stars), 
yield black hole remnants masses to range within 
14-23M© with disk accretion rates in the range: 
0.05 — lM s _1 . If other effects (e.g. disk winds, 
magnetic fields) are important, these masses could 
be lower. But not all stars will form GRBs. For 
stars with insufficient angular momentum to pro- 
duce these disk accretion rates, weak or no explo- 
sions are produced, and the remnant can be much 
more massive (up to the mass of the star). 

This work has been funded by DOE SciDAC 
grant number DE-FC02-01ER41176, a LANL- 
based ASCI grant, and NASA NAG5-9192. 
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Fig. 1. — Critical densities along the pole ver- 
sus effective luminosity [(A:T I y ( ,/4MeV) 2 L l/ ] for set 
values of the annihilation deposition energy (e). 
At high neutrino luminosities, the critical den- 
sity rises sharply as the "Eddington" luminosity 
for neutrinos is reached. The data (e at 20 km 
and L v ) for a series of disks from Popham ct 
al. (1999) are placed on for comparison (tes- 
salated hexagons: M D = 0.01M Q s _1 , triangles: 
M D = O.IMqs" 1 , squares: M D = lM s _1 , cir- 
cles: M D = IOMqs- 1 ). 
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Fig. 2. — Critical polar accretion rates (M) for 
disk accretion rates (Md) of O.IMqs" 1 (dotted 
lines) and lOM s _1 (dashed lines) versus mass 
coordinate of the accreting mass shell for a range 
of black hole rotation values. Note that the ac- 
cretion rate through the disk need not be related 
to (and at late times, is much higher than) the 
accretion along the pole. For comparison, we plot 
the accretion rates of stars ranging from 20-60 Mq 
(Rauscher et al. 2002). 
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Fig. 3. — Accretion rates versus mass for rotat- 
ing and non- rotating 40,60 M Q stars (Heger 2002), 
along with the critical disk accretion rates from 
Fig. 2. 
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